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Resumen 


Esta investigación se desarrolló en la subcuenca del río Puyo, Provincia de 
Pastaza, Ecuador, en los años 2015 y 2018, y tuvo como objetivo 
caracterizar la calidad ambiental de agua en diferentes sectores de la 
subcuenca del río Puyo, incluyendo los principales afluentes. Se clasificaron 
las zonas en alta, urbana y baja. Se realizó un monitoreo en el cauce 
principal y sus afluentes en 28 puntos, registrándose valores de 
conductividad, pH, temperatura, saturación de oxígeno, demanda química 
de oxígeno, demanda bioquímica de oxígeno, y coliformes totales y fecales; 
se aplicaron 413 cuestionarios para identificar los principales tipos de 
actividades que se desarrollan en las zonas aledañas al curso principal y 
afluentes. Los resultados mostraron diferencias en las concentraciones de 
los distintos parámetros en las fechas y puntos de monitoreo, reportando 
mayores concentraciones de DBO», coliformes totales y fecales en las zonas 
del área urbana. El análisis de correlación mostró relaciones significativas 
entre las variables estudiadas, con excepción de la saturación de oxígeno. 
El análisis de varianza mostró diferencias significativas entre las zonas para 
la conductividad eléctrica, y concentración de coliformes fecales y totales. 


Los resultados revelan que la mayoría de aguas residuales se depositan de 


380 
Tecnología y ciencias del agua, ISSN 2007-2422, 12(3), 379-417. DOI: 10.24850/j-tyca-2021-03-10 


Ed 2021, Instituto Mexicano de Tecnología del Agua 
Tecnología y Open Access bajo la licencia CC BY-NC-SA 4.0 
CienciaszAgua 

(https://creativecommons.org/licenses/by-nc-sa/4.0/) 


forma directa en el sistema de alcantarillado y en el cauce principal del río. 
Se concluye que el comportamiento de las diferentes variables registradas 
difiere entre el cuerpo principal y afluentes, incrementando las 
concentraciones de coliformes totales y fecales a partir de la zona urbana 
en el cuerpo principal, y difiriendo el comportamiento del resto de 


parámetros en las tres zonas. 


Palabras clave: contaminación hídrica, agua, población. 


Abstract 


This research was developedin the sub-basin of the Puyo River, Province of 
Pastaza, Ecuador, in the years 2015 and 2018 had as objective: To 
characterize the environmental quality of water, in different sectors of the 
sub-basin of the Puyo River, including the main tributaries. The zones were 
classified as high, urban and low. A monitoring was carried out in the main 
channel and its tributaries in 28 points, registering values of conductivity, 
pH, temperature and oxygen saturation, chemical oxygen demand, 
biochemical oxygen demand, total and fecal coliforms, and 413 surveys 
were applied to identify the Main types of activities that take place in the 
areas surrounding the main course and tributaries. The results showed 
differences in the concentrations of the different parameters in the dates 
and monitoring points, reporting higher concentrations of BOD5, total and 
fecal coliforms, the areas of the urban area. The correlation analysis showed 
significant relationships between the variables studied except for oxygen 
saturation. The analysis of variance showed significant differences between 
the zones for electrical conductivity, concentration of fecal and total 
coliforms. Surveys reveal that the majority of wastewater is deposited 
directly into the sewer system and the main channel of the river. It is 


concluded that the behavior of the different registered variables differs 
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between the main body and tributaries, increasing the concentrations of 
total Coliforms, fecal from the urban area in the main body and differing the 
behavior of the rest of the parameters in the three zones. 


Keywords: Water pollution, water, population. 


Recibido: 25/03/2019 
Aceptado: 25/07/2020 


Introducción y objetivos 


El manejo del agua históricamente se enfocó en suplir las necesidades 
humanas para alimentación y riego, conduciendo a un descuido en la 
disposición de agua fresca para especies y ecosistemas (Ritcher, Mathews, 
Harrison, 8 Wigington, 2003), donde los ríos también son fuente de 
recreación y esparcimiento. La urbanización presenta efectos en 
ecosistemas fluviales en zonas ribereñas, generando cambios en la 
estructura del lecho y aumento de la carga de contaminantes; esto afecta 
las funciones ecosistémicas del curso hídrico. Otras consecuencias son la 
impermeabilización de las áreas de drenaje y captaciones en las cabeceras 
de los ríos (Rórig et a/., 2007). 


El cantón Pastaza, ubicado en la provincia de Pastaza, Ecuador, tiene 


como principal fuente hídrica el río Puyo, el cual presenta influencia de 
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acuerdo con las proyecciones al año 2015 en el cantón Pastazaa 33 550 
habitantes de la zona urbana (Instituto Nacional de Estadísticas y Censos, 
2010); en las zonas rurales, a 5 970 habitantes de la parroquia Tarqui 
(GADPR Tarqui, 2014); 1 027 habitantes de la parroquia Fátima (GADPR 
Fátima, 2015); 232 habitantes de la parroquia Pomona (GADPR Pomona, 
2015); 1 813 habitantes de la parroquia Veracruz (GADPR Veracruz, 2014) 
y del cantón Mera; 10 310 habitantes de la parroquias Shell (GADPR Shell, 
2015), y 1 159 habitantes de la parroquia Madre Tierra (GADPR Madre 
Tierra, 2015); se identifican en toda el área de influencia de la subcuenca 
del río un aproximado de 60 130 habitantes. 


El río presenta dos ramales principales, conocidos como Puyo y Pindo, 
los cuales se unen en el sector Unión Base. En su cuenca media se 
desarrollan actividades turísticas, deportivas y de conservación de 
ecosistemas; atraviesa zonas rurales y urbanas, donde las aguas residuales 
no tienen tratamiento; esto incide en sus características físico-químicas e 
incrementa la importancia de conocer la calidad ambiental de las aguas de 
este curso hídrico, y las relaciones existentes entre las variables físico- 
químicas y microbiológicas con los caudales presentados. 


El presente trabajo se desarrolló en Ecuador, provincia de Pastaza, 
cantones Mera y Pastaza, en la subcuenca del río Puyo, del cual la población 
del cantón Pastaza se abastece para sus requerimientos de agua, donde en 
su curso principal se presentan dos diques ubicados en los ramales Puyo y 
Pindo, y también en sus afluentes donde se localizan cuatro diques 
turísticos. Su cuenca presenta un área de 355 km? y una longitud total de 
cauce de 41.4 km. Además, en varios sectores de su curso principal y 
afluentes se encuentran balnearios naturales donde se realizan actividades 


recreativas. 
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Este trabajo se efectuó entre los años 2015-2018. Se enfocó 
principalmente en el campo temático de calidad del agua y su objetivo fue 
caracterizar la calidad ambiental de agua en diferentes sectores de la 
subcuenca del río Puyo, incluyendo los principales afluentes. 


Materiales y métodos 


Se realizó un monitoreo del año 2015 al año 2018 en 28 puntos: 10 
ubicados en el curso principal y 18 en sus afluentes (Figura 1), cubriendo 
una longitud de cauce de 9.2 km en el ramal Pindo y 17 km en el ramal 
Puyo. Se tomaron muestras compuestas (Jadhav 8 Purohit, 2008), 
registrando valores in situ de conductividad, pH y saturación de oxígeno 
(AWWA Staff, 2010). Se recolectaron muestras para el análisis de demanda 
bioquímica de oxígeno (DBO5) (Standard Methods No 5210B); demanda 
química de oxígeno (HACH 8000) (Hach, 2014); coliformes totales 
(Standard Methods No 9222B) (Rice, Baird, Eaton, € Clesceri, 2012), y 
coliformes fecales (Standard Methods No 9222 D y 92221). En la Tabla 1 se 
reportan coordenadas geográficas, altura, longitud del cauce y gasto base 


para épocas de precipitación normal y baja precipitación. 
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Figura 1. Mapa de ubicación de puntos de monitoreo. 
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Tabla 1. Coordenadas de ubicación de los puntos de muestreo. 


Características 
> e Altura E ' 4 ade Pendiente 
escripción arroquia e cobertura - 
(msnm) al (%) cauca Baja 
el suelo p 
desde la precipitación 
naciente 


Curso principal 


E Altura y Características de cobertura 
Y Descripción del punto Parroquia Gasto base (m”/s) 
(msnm) del suelo 


bial | 
bois | 
facil 
Ela | 
hada 
bici dia 
bel 
bid] 
bad 


x 


180945 | 9824657 Puente sobre el río Indulama 
176956 | 9824201 Balneario del río Indulama* 
Bi 


170644 | 9833900 Puente del río Sandalias BS +P+ZU 
164396 | 9827279 Cabeceras del río Putuimi* 


| 
En 
: 


174838 | 9818474 Bellavista 


164592 | 9830277 Cabeceras río Chingushimi* BS +P+C 


*Puntos también considerados en prueba piloto BS: bosque secundario; P: pastizal, C: 


a. 
a. 
a. 
a. 


cultivo; ZH: zona urbana: ND: dato no disponible. 


Se agruparon los puntos de monitoreo en tres zonas, considerándose 
como zona alta aquella ubicada en la cuenca alta con pendientes > 5 % y 
donde no existen asentamientos poblacionales urbanos; zona urbana, 
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aquella con asentamientos poblacionales superiores a 1 000 personas; y 
zona baja, aquella ubicada aguas abajo posterior a los asentamientos 
poblacionales. También se registró para referencia las precipitaciones 
diarias y mensuales durante el periodo de estudio con el uso de una estación 
meteorológica automática en el punto Estación Biológica Pindo Mirador, en 
la zona alta. Con los resultados obtenidos se aplicó una correlación de 
Pearson (Femández 8: Díaz, 2001) con el uso del programa SPSS (IBM, 
2013), para determinar en función de cada zona la relación existente entre 
los diferentes parámetros evaluados, donde también se incluyeron las 
variables de la precipitación presente el día de monitoreo y la precipitación 
del día anterior. 


Se determinó el estadístico coeficiente de variación desde los puntos 
de monitoreo y se realizó un análisis de varianza con pruebas de 
comparación de Duncan para los diferentes parámetros en función de las 
zonas, para lo cual se utilizó el programa estadístico INFOSTAT (Di-Rienzo 
et al., 2014). Previo al análisis, los valores de porcentaje de saturación de 
oxígeno se transformaron con la función raíz cuadrada (Herrera, 2013). Los 
valores de los diferentes parámetros se compararon con la normativa 
ecuatoriana referida en el Texto Unificado de Legislación Ambiental 
Secundaria (TULSMA), Libro VI Anexo 1 (Ministerio del Ambiente del 
Ecuador, 2015), considerando los criterios de calidad de agua para uso 
recreativo con contacto primario, y criterios de calidad admisibles para la 
preservación de la vida acuática y silvestre en aguas dulces, marinas y de 


estuarios. 


Se aplicó una encuesta en las poblaciones cercanas al curso principal 
y sus afluentes, para identificarlos tipos de actividades productivas que se 
realizan y la disponibilidad de servicios básicos, que fueron distribuidos en 


función de la densidad poblacional de estos sectores; la encuesta se aplicó 
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a un total de 341 pobladores (25 en la zona alta, 212 en la urbana y 104 
en la baja) y a 72 en los afluentes (cinco en la zona alta y 67 en la baja), 
dando un total de 413 cuestionarios, cuyo número supera el valor calculado 
para el tamaño de muestra de poblaciones finitas (Ecuación (1)), que es de 
165, con un nivel de confianza del 99 % (Morales, 2011), considerando una 
población de 60 130 habitantes. 


Tamaño de muestra: 


NZ? 
A (1) 
(N-1)e* + Z*pq 


Donde n es el tamaño de la muestra; N, tamaño del universo (total 
población); Z, nivel de confianza de la estimación, considerando el 99 % de 
confianza (2.576); p, probabilidad de aceptación (0.5); q, probabilidad de 


rechazo (0.5), y e que es el margen de error (10 %). 


Resultados 


La Figura 2 permite observar las características de precipitación 
recolectadas en la estación meteorológica automática ubicada en la estación 
Pindo Mirador, en la zona de estudio en el periodo de muestreo, con 


precipitaciones anuales con una media de 5 005 mm y un promedio mensual 
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de 417 mm, donde sus meses menos lluviosos presentan precipitaciones 
superiores a los 190 mm; se tienen como épocas de mayor precipitación los 
meses de abril de 2015 y enero de 2017; y de menor precipitación 
septiembre de 2015, enero de 2016 y febrero de 2018. Las características 
de precipitación reportadas para los días de muestreo y día anterior fueron 
en 2015 en los muestreos de febrero: 9 y 11.8 mm; marzo: 0 y 2.2 mm; 
abril: 17.6 y 4.4 mm; junio: 0.2 y 4.4 mm. En 2016, en noviembre: 16.4 y 
18 mm; diciembre: 5.6 y O mm. En 2017, en septiembre: 3.7 y 8.2 mm; 
noviembre: 12.2 y 9.8 mm. En 2018, en septiembre: 7.3 y 14.2 mm, y 9.9 
y 9.3 mm, y en octubre: 4.7 y 8.4 mm. 


Figura 2. Precipitaciones mensuales registradas en la estación biológica 
Pindo Mirador, localizada en el área de estudio. Fuente: datos colectados 


de la estación meteorológica automática. 


La conductividad eléctrica a nivel general en todas las fechas (Figura 
3) presenta mayores reportes en el curso principal, en la zona alta, en los 


diques Fátima, Pambay y las Palmas, donde su media sobrepasa los 100 
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uS/cm?, posteriormente disminuye en las zonas bajas de la cuenca; esto 
puede deberse a la presencia de mayores caudales, que generan mayor 
dilución. En los afluentes se reporta mayor conductividad eléctrica en la 
zona alta en el dique Pambay, así como en la zona urbana, en el dique Las 
Palmas; los valores son menores en los otros puntos de monitoreo. El punto 
dique de Shell es el que reporta mayor coeficiente de variación (Tabla 2). 


Cuerpo Principal Afluentes 


- 
4 
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Figura 3. Valores registrados de conductividad. 


Tabla 2. Coeficiente de variación (CV) de las variables registradas. 
Sat. de Colif. Colif. Conduc- 
oxígeno | totales fecales tividad 
E. B. Pindo 
17.8 162.1 148.6 
Dique de 
17.8 162.0 235.3 
Fátima 
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El potencial hidrógeno (Figura 4a) presentó valores entre 6 y 9 en los 


puntos de monitoreo en el curso principal; se compararon los valores con 
los criterios de calidad admisibles para la preservación de la vida acuática 
y silvestre en aguas dulces que establece valores entre 6.5 y 9, y para 
criterios de calidad admisibles para uso recreativo con contacto primario 
con valores entre 6.5 y 8.3 (Ministerio del Ambiente del Ecuador, 2015). La 
mayoría de puntos presentan valores debajo del límite permisible en junio 
de 2015; en abril de 2015 ocurre lo mismo en los puntos Bellavista Bajo, 
Cabecera del río Sandalias 2, puente sobre el río Chorreras, puente sobre 


el río Indulama, balneario del río Indulama, puente Ilayacu y cabeceras del 
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río Chingushimi; en noviembre de 2017 se registró el valor más bajo de 


todos en las cabeceras del río Sandalias 2. 
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Figura 4. Valores registrados de: a) potencial hidrógeno; b) saturación de 


oxígeno. 
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El límite máximo permisible sobre criterios de calidad de agua para 
preservación de la vida acuática y silvestre en aguas dulces no fue 
sobrepasado en los puntos y fechas de muestreo; pero en el criterio de uso 
recreativo con contacto primario, en el cuerpo principal el balnearios dique 
de Fátima sobrepasa el límite en noviembre de 2017 y septiembre de 
2018(b); en los afluentes, los balnearios Dique de Veracruz y balneario del 
río Indulama superan el límite permisible en septiembre de 2018(b). La 
media general muestra los valores más altos en los diques de Fátima y 


Pambay, los cuales se encuentran cercanos al límite. 


Los niveles de saturación de oxígeno (Figura 4b) estuvieron en todos 
los puntos sobre el límite mínimo permisible, que para ambos criterios 
señalados debe sobrepasar el 80 % de saturación, con excepción de las 
fechas de junio de 2015 en los diques Shell, Los Mecánicos, Veracruz; Las 
Palmas, en los puntos Camal Municipal, Cabeceras del río Sandalias 1 y 2, 
Dos Ríos, balneario del río Indulama, cabeceras del río Putuimi y Bellavista; 
también algunos puntos en noviembre de 2016 y octubre de 2018. Los 
puntos donde reiteradamente se reportan valores por debajo del límite 
mínimo son en el cuerpo principal de las cabeceras del río Puyo, el Camal 
Municipal, barrio La Isla, y en los afluentes de los diques Los Mecánicos y 
Veracruz, y en las cabeceras del río Sandalias 2. Los promedios del curso 
principal registran valores superiores al 80 %, disminuyendo en la zona 
urbana e incrementándose en la zona baja, lo cual, al ser el mismo curso 
hídrico, representan una reoxigenación en las zonas bajas. En los afluentes, 
las menores concentraciones se tuvieron en el dique Mecánicos, y en las 
cabeceras del río Sandalias 2. 


La demanda química de oxígeno (Figura 5a) reportó valores más altos 
en los meses de febrero y abril del 2015, principalmente en los sectores de 
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Bellavista Baja, dique de Fátima y barrio La Isla en el curso principal, y en 
los puntos dique de Veracruz, cabeceras del río Sandalias 2, puente sobre 
rio Indulama, balneario sobre el río Indulama y Bellavista en los afluentes; 
y en septiembre de 2018, en el dique Pambay. Los promedios por punto 
muestran mayores valores en el barrio La Isla y Bellavista; baja en el curso 
principal y en las cabeceras del río Sandalias 2, Bellavista y cabeceras del 
río Chingushimi en los afluentes. Este parámetro no está contemplado en 
la normativa legal ecuatoriana, por lo cual no se compara con límites 


permisibles. 
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Figura 5. Valores registrados de concentraciones de: a) demanda química 


de oxígeno; b) demanda bioquímica de oxígeno. 
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En la demanda bioquímica de oxígeno (figura 5b) se consideró el 
límite máximo permisible para criterios de calidad admisible de la DBOs para 
la protección de la vida acuática no impactada, que corresponde a zonas de 
la cuenca alta, donde el límite es 2 mg/l, y con impacto moderado 
correspondiente a zonas urbana y baja, donde el límite máximo es 6 mg//!. 
Los muestreos indicaron valores superiores a los establecido en la norma 
en el curso principal en 2015 en la zona alta en febrero en las cabeceras del 
río Puyo; en el mes de abril, en el punto Pindo; mientras en el mes de junio, 
en la zona alta, en el punto dique de Fátima y en los afluentes en el mes de 
junio, en la zona alta en el dique del ío Pambay; y en la zona baja, en las 
cabeceras del río Putuimi, Nueva Vida y cabeceras del río Chingushimi. En 
los muestreos de noviembre de 2017, septiembre de 2018 (a y b), y octubre 
de 2018, la mayoría de puntos sobrepasa el límite máximo permisible. Los 
valores promedio en el curso principal y afluentes se encuentran dentro de 
lo establecido en los límites permisibles de la normativa ecuatoriana en las 
zonas urbana y baja, donde se considera una condición de impacto 
moderado para la vida acuática; pero no así en la zona alta, donde se 
considera vida acuática no impactada. A su vez, se reportan mayores 
concentraciones en el curso principal en los sectores Rosario Yacu, dique de 


Shell y barrio La Isla, y en los afluentes en el dique Pambay. 


En los valores de coliformes totales (Figura 6a), el criterio de calidad 
de agua para uso recreativo con contacto primario establece un límite 
máximo permisible de 2 000 en su número más probable/100 ml. En lo 
general, todos los puntos sobrepasan el límite máximo permisible, excepto 
las fechas abril de 2015; septiembre de 2018 (a y b); y octubre de 2018, 
donde la mayoría de puntos presenta valores debajo del límite permisible. 


La media general muestra que sólo el punto dique de Los Mecánicos alcanza 
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el límite máximo permisible. En el cuerpo principal se identifica que la zona 
alta presenta menores concentraciones, y éstas se incrementan a partir del 
paso del curso hídrico por la zona urbana; en los afluentes, al sercada uno 
un curso independiente, no se muestra esta tendencia. En la medida de 
coliformes fecales (Figura 6b), el criterio de calidad de agua para uso 
recreativo con contacto primario establece un límite maximo de 200 
NMP/100 ml. El curso principal presenta valores dentro de los límites; lo 
mismo ocurre en la zona alta, los puntos Pindo y dique de Fátima en los 
meses de marzo y junio de 2015, diciembre de 2016 septiembre y 
noviembre de 2017; así como en la zona urbana, en el punto dique de Shell, 
en los meses de marzo y abril. 
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Figura 6. Concentración de: a) coliformes totales; b) coliformes fecales. 


A partir de este punto se incrementa de forma logarítmica la 
concentración de coliformes fecales en la zona urbana y posteriormente 
disminuye en la zona baja, donde se mantiene sobre el límite máximo 
permisible. Los afluentes, en los meses de marzo y junio de 2015, y 
diciembre de 2016, presentan valores debajo del límite máximo permisible. 
Los valores medios muestran a nivel general valores dentro del límite 
permisible en el cuerpo principal en el punto Pindo y en los afluentes en el 


dique “Las Palmas”. El resto de puntos sobrepasa los límites permisibles, 
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teniéndose las mayores concentraciones en el barrio La Isla, el puente de 
Unión Base y Bellavista baja. Los muestreos realizados en 2018 septiembre 
(a y b), y octubre reportaron valores de 1 NMP. 


En la Tabla 2 se muestran los coeficientes de variación para cada 
parámetro en cada punto, donde se puede apreciar a nivel general que las 
variables con mayores coeficientes de variación son coliformes totales y 
fecales, y demanda bioquímica de oxígeno; a su vez, para cada parámetro 
tenemos los puntos con mayor coeficiente de variación en cuanto a 
saturación de oxígeno y coliformes totales en cabeceras del río Sandalias 1; 
coliformes fecales, en el puente de Unión Base; conductividad, en las 
cabeceras del río Putuimi; potencial hidrógeno, en las cabeceras del río 
Puyo; demanda química de oxígeno, en el dique del río Pambay; y demanda 
biológica de oxígeno, en el Puente del río Sandalias. 


El análisis de varianza entre zonas (Tabla 3) reportó diferencias al 
0.05 entre zonas para los parámetros conductividad eléctrica, y coliformes 
totales en el curso principal; mientras que en los afluentes se reportó para 
conductividad eléctrica y coliformes fecales. La prueba de comparación de 
medias de Duncan (Tabla 3) mostró que en la zona baja la conductividad 
eléctrica en el curso principal presenta diferencias con la zona alta y urbana, 
reportando los valores más bajos de las tres zonas, lo cual puede deberse 
a efectos de dilución por aumento de los volúmenes de agua en este sector. 
Este comportamiento también se observa en los afluentes. En el curso 
principal, la medida de coliformes totales en la zona baja presenta 
diferencias con respecto a las otras zonas. En los afluentes, la concentración 
de coliformes fecales reporta mayores valores en la zona baja. En ambos 
casos puede deberse a efectos de acumulación de materiales depositados 
en las aguas residuales colectadas en las zonas alta y urbana. 
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Tabla 3. Análisis de varianza y prueba de comparación de medias de los 


parámetros evaluados entre zonas para el curso principal y afluentes. 


Valor de P 


Parámetro 
Curso principal 


Conductividad 0.0385 
Saturación de oxigeno 0.5834 
0.0742 
Demanda química de oxigeno 0.3074 
Demanda bioquímica de oxigeno 0.5939 
Coliformes totales 


Coliformes fecales 0.0921 


po] 
TI 


Prueba de comparación de Duncan 


Parámetro Coliformes fecales Conductividad 


Media E. E. 
NMP/100 ml Experimental Media uS/cm?2 Experimental 


Alta 386.24 a 225.95 38.67b 
Urbana 390.2 a 291.7 45.57b 
Baja 1026 b 161.25 26.15 a 


Curso principal 


Parámetro Conductividad Coliformes totales 


E. E. 
Media uS/cm2 | experimental | Media NMP/100 ml | experimental 


Alta 60.85ab 10 313.12 a 10 800.55 
Urbana 65.62 b 9 353,55 
j 47.97ab 10 800.55 


a, b: valores con diferencias significativas para p < 0.05. 
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El análisis de correlación (Tabla 4) muestra relaciones positivas 
significativas al 0.01 entre la saturación de oxígeno con la demanda química 
y biológica de oxígeno en la zona baja, y en todas las zonas, entre 
coliformes totales y fecales conductividad y pH; demanda biológica de 
oxígeno con pH y demanda química de oxígeno en todas las zonas. Se 
muestran correlaciones negativas significativas al 0.05 en la zona urbana 
entre DBO con saturación de oxígeno, coliformes totales, y en la zona baja 
entre coliformes totales con pH y demanda biológica de oxígeno. También 
se identifican relaciones significativas positivas al 0.05 entre conductividad 
con demanda química de oxígeno en la zona urbana y con demanda 
biológica de oxígeno en la zona alta. Muestran correlaciones significativas 
en las tres zonas al 0.01 la concentración de coliformes totales con 
coliformes fecales; entre pH y conductividad en la zona alta, y al 0.05 en 
las zonas urbana y baja. Muestra correlaciones significativas 0.05 entre la 
DQO y ph en la zona urbana; entre DBO y DQO al 0.05 en la zona alta, y al 
0.01 en la zona urbana. También en referencia a la precipitación del día de 
muestreo con la precipitación del día anterior, muestra correlaciones 


significativas al 0.01. 


Tabla 4. Correlación entre parámetros. 


Precip. 
] Colif. Colif. mm Precip. D. R 
Zona |Parámetro Conductividad pH DBO |Sat. 02 D. 
totales | fecales muestreo . 
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.036 

-.326 

-.139 

1 -.002 

.267 .021 | -.201 -.201 -.088 
-.023 


-.077 


* La correlación es significativa en el nivel 0.05 (dos colas). 


**lLa correlación es significativa en el nivel 0.01 (dos colas). 


En cuanto a la precipitación del día de muestreo y del día anterior, en 


la zona alta muestra correlaciones significativas con coliformes totales al 


0.05 el día del muestreo y 0.01 el día anterior, y al 0.05 con la conductividad 


y la precipitación en el día anterior en la zona baja; para este caso, las tres 
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zonas reportan el mismo valor, ya que se trabajó solo con una estación 
meteorológica. 


El resultado de las encuestas mostró que en la actividad agrícola en 
el curso principal (Tabla 5) predominan los cultivos de yuca, papa china y 
caña de azúcar; en la zona baja más, el 40 % de los encuestados reporta 
este cultivo; la producción pecuaria en porcentajes inferiores muestra como 
principales al ganado vacuno (22%), porcino y producción de aves (18%), 
y existe predominancia de aserraderos y lubricadores, en especial en la 
zona urbana. En los afluentes, de la actividad agrícola destaca en la zona 
baja el cultivo de papa china, seguido de yuca y caña de azúcar; estos dos 
últimos destacan en la zona alta. En la producción pecuaria, en las zonas 
alta y baja destaca la producción de ganado vacuno, ganado porcino y 
pollos, siendo este último también importante en las zonas bajas, donde se 
reporta baja cantidad de pequeña industria como aserraderos y paneleras. 
Si bien este estudio no contemplo las áreas dedicadas a estas producciones, 
podemos citar a Abril, Ruiz, Alonso, Torres y Cabrera (2016), quienes 
mencionan en un estudio realizado en la provincia de Pastaza, que en su 
mayoría abarca las parroquias ubicadas en el área de estudio, la producción 
de caña de azúcar reporta 1 403.8 ha, 141.3 ha de yuca y plátano, 284.8 
ha de pastos para manejo de ganado vacuno; 17 fincas reportan tener entre 
11 a 50 animales, y en mismo número de fincas de 2 a 10 cabezas de 
ganado vacuno; en ganado porcino, nueve fincas indican tener hasta 10 
animales, y tres fincas reportan tener de 10 a 100 animales. En la 
producción de pollos existen producciones caseras de hasta 100 pollos, con 
un total de 12 fincas; de 101 a 500 pollos en una finca, y producciones de 


más de 2 000 pollos en dos fincas. 


Tabla 5. Actividades que se desarrollan en la zona y servicios sanitarios. 
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En el curso principal en la zona urbana, el alcantarillado llega a un 30 
% de cobertura en la población; en el resto de las zonas es inferiora un 20 
%. La mayoría de zonas poseen servicio de recolección de residuos, pero 
sólo en la zona urbana predomina la recolección diaria. En el curso principal, 
en las zonas alta y urbana predomina el uso de servicio higiénico y pozo 
séptico; mientras en la zona baja se da mayor uso a la letrina; sólo en la 
zona urbana, el 40 % de la población dispone las aguas residuales en el 
sistema de alcantarillado público; en las zonas alta y baja en su mayoría 
son dispuestas en el río o en esteros. En los afluentes, en la zona alta existe 
mayor uso de servicio higiénico y las aguas residuales se depositan en su 
mayoría al alcantarillado público; en las zonas bajas se usa 
mayoritariamente letrinas y pozos sépticos, y las aguas son depositadas en 
el río y esteros. Los resultados de los servicios básicos con los que cuentan 
las viviendas al ser relacionados con las concentraciones de coliformes 
fecales indican que en las zonas alta y baja que presentan características 
rurales, la deposición de materia fecal presentaría mayor cantidad de 
descargas con poco volumen, generando una dispersión del contaminante; 
mientras en la zona urbana, al ser depositado en el alcantarillado, concentra 
la materia fecal en el sector, recibiendo también aportes de descargas 
dispuestas directamente al río. Esto ratifica que la influencia del tipo de 
zona incide en la concentración de coliformes fecales, siendo más evidente 
en la zona urbana, donde se identificaron diferencias estadísticas 


evidenciadas en el análisis de varianza. 


Discusión 
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El estudio muestra una interacción entre el tipo de zona, donde existe una 
influencia de la población sobre la concentración de los diferentes 
parámetros; donde en la zona alta, al tener una cobertura principalmente 
de bosque secundario y menor cantidad de población, muestra en sus 
parámetros valores más ajustados a los criterios establecidos en la 
normativa ecuatoriana. La zona, al presentar una topografía ondulada, 
muestra mayor cantidad de pendientes, presenta cambios en el relieve en 
a partirde los 8 km desde la naciente en ambos ramales, donde de acuerdo 
con los criterios de Martínez (2006) presentan una pendiente alta (> 7 %) 
y media (de 2 a 7 %), lo cual influye en tener una mayor velocidad en el 
flujo y más presencia de saltos de agua, generando mayor oxigenación; sin 
embargo, el tipo de origen de cada uno de los ramales también muestra 
una respuesta diferente, pues el ramal Pindo nace de la acumulación de 
aguas de un pantano y el punto de monitoreo se encuentra a 4 km de la 
naciente; mientras que el punto cabeceras del río Puyo se encuentra en la 
misma naciente que proviene de fuentes subterráneas. Así, en el primer 
caso tiene un mayor recorrido e influye en tener una mayor saturación de 
oxígeno; sin embargo, en el ramal Puyo, a 8 km de la naciente, el punto 
dique de Fátima muestra valores promedios similares a los de la estación 
Pindo Mirador, lo cual, al considerar que aguas arriba de estos puntos no 
existen asentamientos poblacionales, muestra que distancia y pendiente 
generan estos niveles de saturación con medias del 90 %. Esto también se 
refleja en menores concentraciones de DQO, DBO5 y coliformes totales, 


pero no así en las concentraciones de coliformes fecales. 


Al analizar en el cuerpo principal las concentraciones de coliformes 
totales y fecales en ambos ramales muestran un efecto de acumulación a 
partir de que la corriente atraviesa la zona urbana; al considerar los últimos 
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puntos en esta zona antes de la zona baja-rural, los caudales alcanzan los 
7.42 m3/s en el ramal Puyo y 6.03 en el ramal Pindo; esto nos plantea que 
tales volúmenes de agua se ven superados en su capacidad de dilución de 
las dos variables. 


Romero, Chaile y Amoroso (2010), en estudios en ríos de montaña, 
presentan un incremento en el valor de potencial hidrógeno en las zonas de 
cuenca baja con respecto a las de cuenca alta, donde no se identifica la 
presencia de grandes poblaciones, lo cual no es observado en este estudio, 
e indica que la presencia de población influye en el comportamiento de la 
calidad de agua. 


Pérez y Rodríguez (2008) señalan que cambios en el pH pueden 
indicar el ingreso de fertilizantes, de procesos de eutrofización; Marques, 
Albertoni y Palma-Da-Silva (2012) indican valores altos antes de la 
presencia de represa; en comparación con los resultados del presente 
estudio, en la zona alta en el curso principal reportan valores más elevados 
el dique de Fátima, y en los afluentes en el dique Pambay y en la zona 
urbana en el dique Las Palmas, los cuales tienen características de no 


presentar poblaciones en su curso superior. 


Los valores de oxígeno disuelto son superiores al planteado por Pérez 
y Rodríguez (2008), quienes indican que en aguas tropicales, niveles 
inferiores al 55 % afectan la vida acuática. Esto puede deberse a ser un río 
proveniente de montaña, donde la presencia de mayor cantidad de rápidos 
faculta el intercambio de oxígeno con el medio, aumentando su 


concentración. 


En la demanda biológica de oxígeno, Abril, Rodríguez, Landázuri y 
Santi (2014), en estudios realizados entre marzo y mayo del 2012 en el 


dique de Shell, muestran valores entre 2 y 16 mg/l en este punto, con una 
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media de 9.3 mg/l ml, y muestra que, en las condiciones del año 2015, la 
presencia de sistemas de tratamiento influye en mitigarlas concentraciones 
de dicho parámetro en este punto; en ese periodo, los valores se reportan 
dentro de los límites permisibles. 


Abril et al. (2014), en estudios realizados entre marzo y mayo del 
2012 en el dique de Shell, indican valores de coliformes totales entre 800 y 
1 000 000 de ufc/100 ml, con una media de 257 000 ufc/100 ml, y muestran 
que, en condiciones del periodo evaluado, la presencia de sistemas de 
tratamiento influye en disminuir las altas concentraciones de este 
parámetro en tal punto. A partir de este punto, en lo general, el curso 
principal excede los límites permisibles en las zonas urbana y baja, pero en 
los muestreos que tuvieron mayor precipitación (marzo y abril 2015) 
muestra una tendencia a disminuiren las zonas bajas. Los mismos autores 
señalan que coliformes fecales muestran valores entre 110 y 900 000 de 
ufc/100 ml, con una media de 228 078 ufc/100 ml, e indican que, en las 
condiciones del año 2015, la presencia de sistemas de tratamiento influye 
de manera significativa en reducir las altas concentraciones de dicho 
parámetro en tal punto. Esto también muestra que los sistemas de 
tratamiento eliminan una importante carga de coliformes fecales, y su 


remoción es mayor en proporción a la obtenida en coliformes totales. 


Marques et al. (2012) indican muestreos en zonas previas y 
posteriores al represamiento de agua en el río Antas en Brasil, donde la 
concentración de las demandas biológica y química de oxígeno disminuyen 
en zonas posteriores al represamiento, mientras que las concentraciones de 
coliformes totales y fecales se incrementan; esto se puede observar entre 
los puntos cercanos dique de Shell y Camal Municipal. 


Torres, González Rustrián y Houbron (2013), en monitoreos 
realizados en el río Seco, México, con una extensión del curso de 70 km y 
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una área de cuenca de 438 km?, con precipitaciones en toda la época 
lluviosa de 1 426 mm y en épocas secas de 171 mm, que presenta 
características en su parte alta de pendientes del 15 % hasta el primer 
punto de muestreo y del 3 % en el resto de la cuenca, muestran las menores 
concentraciones de demanda química de oxígeno en el primer punto de 
monitoreo en la parte más alta de la cuenca; lo mismo se puede observar 
en el presente estudio, donde el área de cuenca es de 350 km? y la 
extensión del curso de 40 km, a pesar de presentar mayores 
precipitaciones. 


En la zona baja, el uso de letrinas y la presencia de ganado vacuno, 
porcino y planteles avícolas son factores que aportan cargas de materia 
fecal al río. Es necesario mencionar que la presencia de estas bacterias no 
necesariamente representa la existencia de contaminación fecal humana, 


pues éstos también se encuentran en el suelo y la vegetación (Romero et 
al., 2010). 


Se observa que las concentraciones de DBOs (Whipple, Hunter, € Yu, 
1974) y coliformes fecales (Johnston, Detenbeck, 8 Gerald, 1990) son 
mayores en las zonas urbanas, donde existe mayor población, lo cual se 


vincula con lo citado por Rónrig et al. (2007). 


En la correlación de variables con la precipitación del día de muestreo, 
sólo la concentración de coliformes totales presentó correlaciones 
significativas al 0.01 en la parte alta. Si aplicamos el criterio hidrológico del 
tiempo de concentración, el cual expresa el tiempo de recorrido de un flujo 
de agua entre dos puntos en una cuenca, a través de la fórmula de Kirpich 
(Gaspari, Rodríguez, Senisterra, Delgado, 8, Besteiro, 2013), el punto dique 
de Fátima presenta un tiempo de 1.9 horas y el de la Estación Biológica 
Pindo Mirador de una hora, lo cual fundamenta que estos puntos ubicados 


en la parte alta presenten influencia y correlaciones significativas con 
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respecto a la precipitación del día de muestreo, al recibir aportes de 
caudales generados por precipitaciones en horas previas. 


Conclusiones 


La topografía, el uso de suelo y la presencia de poblaciones influyen en las 
características de las concentraciones de los parámetros saturación de 
oxígeno, DQO, DBOD5, coliformes totales y fecales en el cuerpo principal y 
afluentes. 


En el cuerpo principal, las concentraciones de coliformes totales y 
fecales se incrementa en escala exponencial a partir de la zona urbana; esto 
se refleja en sus promedios, sobrepasando los 10 000 NMP/100 ml en 
ambas medidas; al tener mayores aportes de aguas residuales dispuestas 
en el sistema de alcantarillado y posteriormente vertidos en el curso hídrico, 
sobrepasa el límite establecido en la normativa. Esta tendencia no se 


observa en los afluentes. 


Los coeficientes de variación indican que las variables concentración 
de coliformes totales y fecales reportan mayor vaniabilidad, con rangos de 
100.6 a 274 en coliformes totales, y de 88.1 a 250.9 en coliformes fecales 
para los puntos de monitoreo; a su vez, el pH es el que menores coeficientes 
de variación reporta. El análisis de correlación muestra al DBO con mayor 
cantidad de correlaciones con el resto de variables, principalmente con la 


DQO, donde en la zona alta y media presenta correlaciones significativas al 
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0.01; a su vez, la precipitación del día de muestreo significativa al 0.05 y 
día anterior significativa al 0.01, presenta relación con la concentración de 


coliformes totales en la zona alta. 


El análisis de varianza mostró diferencias para p < 0.05 en la 
conductividad eléctrica para el curso principal y afluentes, teniéndose 
menores valores en la zona baja, con una media de 47.97 uS/cm? en el 
curso principal y 26.15 uS/cm? en los afluentes; al estar esta variable 
relacionada con las concentraciones de sales disueltas, muestra un efecto 
de dilución para este parámetro. Las concentraciones de coliformes fecales 
en los afluentes y totales en el curso principal presentan los mayores valores 
en la zona baja. 


La población y el tipo de disposición que dan a sus aguas residuales 
son el principal aporte que influye en las concentraciones de DBO, DQO, y 
coliformes totales y fecales, pues en las zonas urbana y baja más del 70 % 
lo disponen al alcantarillado y río, y sólo en los sectores aledaños al dique 
de Shell se tienen plantas de tratamiento, lo cual también se observa en el 
análisisde varianza para esta variable, lo que genera un efecto acumulativo 
hacia la zona baja, donde se reporta mayor concentración, con una media 
de 48 833.3 NMP/100 ml. 


Las labores agrícolas y pecuarias se desarrollan principalmente en la 
zona baja-rural, donde en los afluentes, con base en al análisis de varianza, 
se reportan mayores concentraciones de coliformes fecales, con una media 
de 1 026 NMP/100 mil. 


A nivel general se recomienda que en estrategias para el manejo del 
recurso hídrico en zonas similares a la del estudio, los sectores ubicados en 
la parte alta antes de las poblaciones se deben enfocar en la preservación 


de los bosques; en la zona urbana, el enfoque debe direccionarse al 
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tratamiento de las aguas residuales domésticas; y en la zona baja, hacia 
las actividades agrícolas y las aguas residuales domésticas. 
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